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Abstract. The article begins with a quantitative version of the martingale 
central limit theorem, in terms of the Kantorovich distance. This result is then 
used in the study of the homogenization of discrete parabolic equations with 
random i.i.d. coefficients. For smooth initial condition, the rescaled solution 



(N 

o 

*H ' of such an equation, once averaged over the randomness, is shown to converge 

polynomially fast to the solution of the homogenized equation, with an explicit 
exponent depending only on the dimension. Polynomial rate of homogenization 
for the averaged heat kernel, with an explicit exponent, is then derived. Similar 
l/j ' results for elliptic equations are also presented. 
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1. Introduction 



1.1. Main results. The main goal of this article is to give quantitative estimates in 
the homogenization of discrete divergence form operators with random coefficients. 
Writing B for the set of edges of Z d , we let u = (w e ) e gB be a family of i.i.d. random 
variables, assumed to be uniformly bounded away from and infinity, and whose 
t^J- ' joint distribution will be written P (with associated expectation E). The operator 

whose homogenization properties we wish to investigate is 

(i.l) #7(*) = X>x,v(/(v) - /(*)) (xez% 

where we write y ~ x when x, y G Z are nearest neighbours. For a bounded 
continuous / : R d — >• R, we consider u^ the solution of 

^ ' ( fl»#(e) 

b ' ^— = L u u& on R+ x Z d , 

(DPE^) <^ dt 

u (e) (0,-) =/(e •) onZ d , 

and u e (t, x) = u^ 6 ' (e^ 2 t, [e ^ 1 x\ ). There exists a symmetric positive-definite matrix 
A (independent of /) such that the function u £ converges, as e tends to 0, to the 
function u solution of 



(CPE) 



— = -V • (AVu) on (0, +oo) x M d , 




The notions of being a solution to (DPE" ) or (|CPE[) . and of the convergence of 



u E 



to u, will be made precise later on. For a bounded measurable function / : M. d 



l 
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and a = {ct\, ■ ■ ■ , ay) & N rf , we call 

d xf 1 ...x" d f 



d\<*Wf 



(Mi = E 7 "j) 



dx" 1 • • • dx a d d 

a weak derivative of order |a|i, where the derivative is understood in the sense of 
distributions. 

Here and below, we write |_a;J for the integer part of x, a Ab = min(a, 6), a V b = 
max(a, b), log, (x) = log(x) V 1, and |£| for the I? norm of £ G W l . The main 
purpose of this paper is to prove the following theorems. 

Theorem 1.1. Let m = [d/2\ + 3 and S > 0. There exist constants C$ (which 
may depend on the dimension) and q such that, if the weak derivatives of order m 
of f are in L 2 (M. d ), then for any e > 0, t > and x £ R d , one has 

(1.2) \E[u 6 {t,x)]-u{t,x)\ 

d I d 

^R./ius+ca ( t + ^) ii/ii2+Eii%'/n 



^1,5 I - 



where 
(1.3) 



*,,j(«) 



m 1 / 4 i/d = l, 

log^u" 1 ) u 1 / 4 ifd = 2, 
u l ' 2 ~ s ifd^i. 



Theorem 1.2. Letp^(x,y) be the heat kernel associated to L' 

1 



Pt&,y) 



(2irt) d / 2 VdetA 



exp 



— (y~x) T A (y-x) 



be the heat kernel associated to -k V • AW , and S > 0. There exist constants c > 
(independent of 5), q, Cg, eg > swc/i that for any e > 0, t > satisfying ej \ft sj e^ 
and any a; (E R d , one /las 



(1.4) 



,-d 



E[p-_ 2i (0, Le" 1 ^)] -P,(0,^)| 



C<5 / T / 



l/(d+3) 



exp 



(^-A|e x x| 



In particular, for any s > 0, one has 

sup sup|e~ d E [p%- H % [e-^J)] -ft(0,s)| - O ((*^(e 2 )) V(d+3) ) 



as e tends to 



Remark 1.3. For a given smooth function / and a fixed i > 0, the r.h.s. of (|1.2|) is 
of the order of 

y/e if d = 1, 

log^e" 1 )^ ifd = 2, 



(1.5) 



e 1 - 5 ' if d Si 3, 



where <5' = 2<5 > is arbitrary. Similarly, for fixed £ and x, the r.h.s. of (|1.4p is of 
the order of 

a/8 



(1.6) 



:' ifd=l, 

log^e" 1 ) e 1 / 10 ifd = 2, 
£ V(d+3)-5" if O 3, 



where 5" = 26/ (d + 3) > is arbitrary. 
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1.2. Context. Homogenization problems have a very long story, going back at 
least to |Mal ISt] . Rigorous proofs of homogenization for periodic environments were 
obtained in the 70's (see for instance [BLP] or [JKOI Chapter 1] for references), 
and for random environments with [Ko78, Yu80, PV81, Kii83 . Classical methods 
used to show homogenization typically rely on a compactness argument, or on the 
ergodic theorem, both approaches leaving the question of the rate of convergence 
untouched. 

For continuous space and periodic coefficients, [JKOI Corollary 2.7] uses spectral 
methods to show that 



-<! 



p^ H (0 1 e- 1 x)-p t (0,x)\^C 



t (d+l)/2- 



For random coefficients, available results are much less precise. For continuous 
space, |Yu86j gives an algebraic speed of convergence of u e to u for the elliptic 
problem and d ^ 3, without providing an explicit exponent. In |CS10j . the much 
more general case of fully nonlinear elliptic equations is considered, and a speed of 
convergence of a logarithmic type is proved. 

Here, we focus on the convergence of the average of u £ to u. This approach 
has been considered in [CSllj for the elliptic problem. There, it is shown that 
the suitably rescaled Green function, once averaged over the randomness of the 
coefficients, differs from the Green function of the homogenized equation by a power 
of e. The exponent obtained is implicit, and depends on the ellipticity condition 
assumed on the random coefficients. 

In contrast, Theorems 11.11 and 11.21 provide explicit exponents, that depend only 
on the dimension. I conjecture that the correct order of decay with e in Theorem ll.il 
should be 

y/e if d = 1, 

log^ 1 ) e ifd=2, 

e if d ^ 3. 

This differs notably from what is obtained in Theorem 11.11 only when d = 2. On 
the other hand, it may well be that the assumption of high regularity on / is only 
an artefact of the methods employed. 
The fact that 

(1.7) Bupsup|e- ,, i£_ at (0,Le" 1 3j)-p t (0,s)| -2±+ 

is known at least since JBH09] , where the much more difficult case where the random 
coefficients are Bernoulli random variables is considered (in this context, the heat 
kernel should be considered only within the unique infinite percolation cluster). 
Yet, for strictly positive random coefficients, this convergence does not hold if the 
distribution of the random coefficients is allowed to have a fat tail close to and 
when d ^ 4 [BBHK08] IBB10) . Under the same circumstance and when p" is 
replaced by its average in (11.7[) . the convergence fails to hold in any dimension 
FM06] (see however (5.BDH10, Proposition 7.2] for a nice way to get around this 
problem) . 

An evaluation of the gap between the average of u e and u naturally calls for 
estimates on the size of the random fluctuations of u e around its average. In 
this direction and for the elliptic problem, [CNOOaj obtains algebraic decay of the 
variance of u £ (integrated over space). The exponent obtained is implicit, and 
depends on the ellipticity condition^. 



A. Gloria has announced improved estimates on this problem. 
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1.3. Our approach. In order to prove Thcorcm ll.il we will use the representation 
of u e as the expected value over the paths of a random walk, that we write (X t )t^o- 
This random walk has inhomogeneous jump rates given by the (w e )eeB, and L u is 
its infinitesimal generator. For instance, one has 

u e (t,0)=E%[f(eX e - H )}, 

where we write Pq for the distribution of the random walk starting from 0, and 
Eq for its associated expectation. From the PDE's perspective, this approach can 
be seen as the method of characteristics, except that here those characteristics are 
random. 

From the probabilist's perspective, the (pointwise) convergence of u e to u is 
equivalent to the claim that the random walk, after diffusive scaling, satisfies a 
central limit theorem. Quantitative estimates should thus follow if one can provide 
with rates of convergence in this central limit theorem. 

In [Mollj , it is shown that there exist constants C, q ^ such that for any £ of 
unit norm, 



(1.8) 



sup 



Z-x t 



€ x 



$(x) 



^C 



r i/io 




if d = 


= 1 


log 9 + (i) t- 


-1/10 


if d = 


-2 


log+(t) f 


-1/5 


if d = 


:3 


r i/s 




if d> 4 



where $ is the cumulative distribution function of the standard Gaussian random 
variable, and cr(£) = £ • A£. 

This result has two important weak points: (1) the rates are far from the usual 
t^ 1 ' 2 one obtains for sums of i.i.d. random variables, and (2) the theorem only gives 
information about the projections of X t onto a fixed vector. We shall find ways to 
overcome these two problems. 

The classical approach for the proof of a central limit theorem for the random 
walk consists in decomposing it as the sum of a martingale plus a remainder term, 
and then show that the martingale converges (after scaling) to a Gaussian random 
variable, while the remainder term becomes negligible in the limit. 

In view of this, what should be done is clear: we should first find a quantitative 
estimate on how small the remainder term is, and second, show that the martingale 
converges rapidly to a Gaussian. This is indeed the method used in [Moll . The 
control of the remainder term given there is satisfactory, and the problem lies with 
the quantitative central limit theorem for the martingale part. 

This quantitative central limit theorem relies on the fact that one can have a 
sharp control of the variance of the quadratic variation of the martingale. It is 
shown that, after scaling, this variance decays like £ _1 when d ^ 4, which is the 
best possible rate. However, given such a control, the quantitative CLT (due to 
[HB701 IHa88j ) used there only yields a decay of £ -1 / 5 in this case. 

Surprisingly, this exponent 1/5 is best possible |Mol2j . To overcome this ob- 
struction, we derive new quantitative CLT's for martingales, that will not yield a 
Bcrry-Esseen type of estimate, but rather measure 

i-x t 



sup 



EEn 



/ 



r(0Vi) 



/d4> 



where £ is a class of functions. When £ is the class of bounded 1-Lipschitz func- 
tions, the supremum is often called the Kantorovich(- Rubinstein) distance. We also 
consider £ to be the class of bounded C 2 functions that have first derivative bounded 
by 1 and second derivative bounded by k, and call it the k-Kantorovich distance. 
The martingale CLT's obtained hold for general square- integrable martingales, and 
are of independent interest. 
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Once equipped with these martingale CLT's, we apply them to the one-dimen- 
sional projections of the random walk (X t ), and for d ^ 3, we obtain rates ap- 
proaching the i.i.d. rate of f 1 !" 1 . To do so, we use estimates derived in |Mollj . 
most importantly on the variance of the quadratic variation of the martingale. 
These in turn are consequences of the LP boundedness of the corrector (for d ^ 3, 
and with logarithmic corrections for d = 2), and of a spatial decorrelation property 
of this corrector, proved in |GQ111 Theorem 2.1 and Proposition 2.1]. 

In order to obtain Theorem ll.il we need to carry the information obtained on the 
projections of X t to X t itself, in a kind of quantitative version of the Cramer- Wold 
theorem. This is achieved through Fourier analysis, at the price of requiring the 
existence of weak derivatives of higher order. 

The key observation that enables to go from Theorem 11.11 to Theorem 11.21 is the 
high regularity of the averaged heat kernel. In contrast to the true heat kernel, the 
averaged one has a gradient which is bounded by a constant times the gradient of 
p, as is proved in jCNOOal IHD05] . 

The estimates due to |G011] are the only place where the assumptions of inde- 
pendence and uniform ellipticity of the coefficients come into play. In particular, if 
it is shown that these estimates are valid for certain correlated environments, then 
the present results automatically extend to this context. The present results should 
also extend to continuous space with only minor change, as long as the estimates 
of [GOllj remain true in this setting^. 

1.4. Organization of the paper. We introduce the (fc-)Kantorovich and Kol- 
mogorov distances in section [5] In section [31 we consider general square- integrable 
martingales, and derive quantitative CLT's with respect to the (fc-)Kantorovich 
distances. We then apply these results to projections of the random walk X t in 
section 0J The homogenization setting is taken up in section [SJ and Theorem 11.11 
is proved. Theorem 11.21 is then derived in section [BJ Finally, similar results for the 
homogenization of elliptic equations are presented in section [7) 

2. Distances between probability measures 

A function / : R m — > 1" is said to be fc-Lipschitz if for any x,y £ R m , one has 
\f(y) ~ f( x )\ ^ k\y — x\. Let v, v' be probability measures on R, and let F Vl F v i 
be their respective cumulative distribution functions. We define the Kantorovich 
distance between v and v' as 



(2.1) di(z;,i/) =sup 



Jfdu-Jfdu' 



, / bounded and 1-Lip. > , 



and the Kolmogorov distance between v and v' as 

(2.2) (!«,(!/, i/) = sup \F v ,{x) - F v {x)\ = || F„, - F u \\ x . 

The notation for the Kantorovich distance becomes more transparent once we notice 
that (see for instance |Vil Theorem 1.14 and (2.48)]) 

(2.3) di(i/y)= [\F„,(x)-F v (x)\dx = \\F v ,-F v \\ 1 . 



As we will see below, bounds in the martingale CLT are improved when measured 
with the Kantorovich distance instead of the Kolmogorov distance. We now intro- 
duce weaker forms of the Kantorovich distance, for which the rates of convergence 



A. Gloria and F. Otto have indeed announced extensions of their results to the continuous 
setting, and to random environments satisfying a Dobrushin-Shlosman mixing condition. 
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will be even better. For any k G [0, +00], we define the k-Kantorovich distance as 



di,fe(^,^') = sup 



/di/ - / fdi/ 



, /ec?(R,R),||/ / || 00 <i J ||/ ,, || 00 <fcL 



where C^ ( 



L) is the set of bounded twice continuously differentiable functions from 
R to R. For fc < fc', one has di, fe < di,*/ s£ di !OC = di. Note that if / e C^(R,R), 
then 



(2.4) 



fdv - / /d!/ 



<ll/'l|oodi, 



3 /ll/' 



>y). 



In the sequel, if X follows the distribution v and y the distribution 1/ , we may 
write di(X, V) to denote di(^, 2/), or also di(X, F„<) if convenient. If X and Y are 
defined on the same probability space with probability measure P and associated 
expectation E, then for any 1-Lipschitz function /, we have 

\E[f{X)\ - E[f(Y)]\ ^ E\f(X) - f(Y)\ < E\X - Y\, 

and hence 

(2.5) d x (X,Y)^E\X-Y\. 

Similarly, if X follows the distribution v and Y the distribution z/, we write 
di } h(X,Y), d\ t k{X,F v t) or di,*;^, v') as convenient. 

3. Martingale CLT 

For a square-integrable (cadlag) martingale (^t)te[o,il defined with respect to 
the probability measure P and the (right-continuous) filtration {Ft)t~^o, we write 
((M)t)te[o,i] for its predictable quadratic variation, 

AM(t) = M t - lim M s , 



and 



L 



■2l> 



E 



E |AM(t)| 2 ^ 



Recall that we denote by $ the cumulative distribution function of the standard 
Gaussian random variable. In |Ha88] , the following is proved. 

Theorem 3.1 ( Ha88 ). For any p > 1, there exists C p (independent of M) such 
that 



(3.1) 



cUMx, $) ^ c p (l^ (2p+1) + MM), - 1\\;/^) . 



Our first result consists in showing that one can get sharper bounds if one replaces 
the Kolmogorov distance by the (fc-)Kantorovich distance in fl3.1[) . 

Theorem 3.2. For any p > 1, there exists C p (independent of M) such that 

(3.2) d!(M 1} *) < C p L 1 2 / p {2p+1) + 2\\(M) 1 - l\\l /2 , 
and for any k Js 0, 

(3.3) d 1)fc (M 1} 3) < C p ^ (2p+1) + ^ p + (fc V 1)11(71/)! - l||x. 

Remark 3.3. Naturally, one has ||(M)i — l||i ^ ||(M)i — l|| p , and the statements are 
only interesting when this quantity, and also L2 P , are small, so Theorem l3.2l indeed 
provides better rates of convergence than Theorem l3.il It is shown in [Mo 12] that it 
is not possible to change the exponent p/[2p+\) appearing on the term ||(M)i — l|| p 
in the r.h.s. of (|3.1|) by any higher exponent. It would be interesting to investigate 
how sharp (I3.2p is in this respect. The term ||(M)i — l]|i appearing on the r.h.s. 
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of (|3.3[) cannot be improved. Indeed, let (B s ) s ^,q be a standard Brownian motion, 
and consider the martingale M s = Bn+ S \ s . Since the martingale is continuous, L2 P 
vanishes, while one has ||(M)i — l||i = s. On the other hand, the cosine function 
has first and second derivatives bounded by 1, and thus 

d M (Mi, *) > S[cob(Bi)] - £[cos(Afi)] = e" 1/2 - e -( 1+£ )/ 2 ~ £ ( e -> 0), 

thus justifying optimality of the exponent on ||(M)i — l||i. 

Remark 3.4. A quantitative martingale CLT expressed in terms of the Kantorovich 
distance was already formulated in [RRj Theorem 8.1.16]. Terms involved in the 
bound are however difficult to estimate in practical situations, in contrast to what 
is obtained in Theorem [ 



In order to prove Theorem l3.2l we will rely on the following non-uniform version 
of Theorem O 

Theorem 3.5 ( HJ88, Ha88]). For any p > 1, there exists C p (independent of M) 
such that if Li p + \\ {M)\ — 1||^ ^ 1, then for any x E M. 

\P[ Ml < x] - $(,)! < ^L- (l%^ + MM), l||f 2 ^>) . 

HJ88] Theorem 1] is the equivalent statement concerning discrete time mar- 
tingales. Theorem 13.51 can be derived from its discrete time version by apply- 
ing the approximation procedure explained in |Ha881 Section 4] (in [Ha§8], lo- 
cally square-integrable martingales are considered, while we stick here with plainly 
square-integrable martingales. There is no loss of generality however, since a lo- 
cally square-integrable martingale is in fact a square-integrable one if ||(M)i — l||i 
is to be finite. One can thus skip the localization procedure at the end of |Ha881 
Section 4]). 

Proof of Theorem \3.2\ We start by proving that there exists C p (independent of 
M) such that (J3.2I) holds. We decompose the proof of this into three steps. 

Step 1.1. We first prove the claim assuming that 

(3.4) (M) 1 = 1 a.s. 

and that L>2 P ^ 1. Under this condition, Theorem 13.51 ensures that 

|P[M 1 ^]-$(z)|< I ^L*f^. 

We thus have, after possibly enlarging C p , 

di(Mi, $) = J \P[M X ^x]- §(x)| dx < C P L\ 

which is the desired result. 

Step 1.2. We now no longer impose that condition (J3.4D holds, but keep with the 
assumption that Li p $J 1. Following an idea probably due to [Dv72 , we introduce 

t = sup{£ < 1 : (M) t < 1}. 

Note that r is a stopping time, since 

{r^t}=f]{(M) t+e >l}eT t . 

We define 

(M) 1 if(M) 1 <l, 

lim t _j. r - (M)t otherwise, 



-l/(2p+l) 
J 2 P 



(M) T - 
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M T - = 



and 

Mi if (M) i^l, 

lim 4 _j. T - M t otherwise. 

Note that (M) T - ^ 1. Let (B s ) s ^q be a standard Brownian motion, independent 
of the martingale. We define 



M,= 



M s 
M T - 

M T -+B s _t 
M T - + B x _ {m)t _ 

By construction, M is a martingale, and 

(M) 2 -(M)i 



if o «; s < t, 

if r <s < 1, 

if l<a<2- (M) T -, 

if 2 - (M) T - < s s: 2. 



1 - (M) 7 



hence (M) 2 = 1. Naturally, the fact that M is defined on [0,2] instead of [0,1] 
plays no role, and this martingale satisfies condition (J3.4I) (at time 2). Let us write 



■>2p 



E 



E i aa? wi 



2l> 



0<t<2 



< L 2p < 1. 



We learn from the the first step of the proof that 

(3.5) dl (M 2 , $) < C p L\> p ^ < (7 p L 2 f^. 
We now want to use the fact that 

(3.6) di(Mi,$) ^di(Mi,M 2 ) + di(M 2 ,$) 
to estimate di(Mi, $). In view of (|2.5[) . we have 



di(Mi,M 2 ) s^E 



Mi -M 2 



Note that 



and thus 



Mi - M 2 = M x - M r + AM(r)l (M>I> i - £i_<m> t _ , 



Mi -M 2 



< £ [|Mi - M r |] + E [|AM(r)|] + E 



B l-(M) T . 



Let us write a\ +a 2 + 03 for the latter sum, with obvious identifications. We bound 
the contribution of each of these terms successively. 

ai < E [(Ah - M T ) 2 ] V2 = E [<M)i - (A/) r ] 1/2 , 

since r ^ 1 is a stopping time. Now, either r = 1, in which case (M)i — (M) T = 0, 
or t < 1, in which case (M) r ^ 1. In both cases, we have 

(M)i - <M) T < |<M)i - 1| , 

and thus 

ai^lKMh-lH}/ 2 . 

As for a 2 , we have 



(3.7) a 2 = £ [|AM(r)|] < £ [|AM(r)| 2p 

For the third term, we have 



V(2p) 



^i 



l/(2p) 
2p 



as = c E [|1 - (M) r -| 1/2 ] < c £ [|1 - (Af) r -|] 1/2 , 
where c = E'fl-Bil] ^ 1. We decompose the last expectation as 

E[\l-(M) T -\ l {M)l<1 ]+E[\l-{M)r-\ l(M)x>l] 
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The first term is bounded by ||1 — (M)i||i, while the second term is smaller than 

£[A(M)(T)]a[(AM(r)) 2 ] 

(to see this, consult for instance the proof of [JS1 Theorem 4.2]). The latter is 
bounded by 

B[(AM(r)ff /p aJ. 
To sum up, we have shown that 

6 1 (M 1 ,M 2 ) < IKA/h - 1||} /2 + L\' p {2p) + j\\l- (MjilU + L%* 
(3.8) < 2\\{M) l -l\\\ /2 + 2L\' p [2p) . 

Since we assume that L 2p ^ 1, we have L 2 L < L 2p p , and equations (|3.6[l . 
(|3~5)1 and (O give us that 

diCMx,*) < (<5„ + 2)L 2 f p+1) + 2||<M)i - 1||} /2 , 
which is what we wanted to prove. 

Step 1.3. There remains to consider the case when L 2p > 1. It follows from (|2.5[) 
that 

d 1 (M 1 ,^)^c+\\M 1 \\ 1 , 
where c is the L\ norm of a standard Gaussian, c ^ 1. Moreover, 

||Mi||i < HMxHa = \\(M)i\\l /2 < (1 + ||{M)i - lh) 1/2 . 

As a consequence, it is always true that 

di(M 1) $)<2+||(M) 1 -l||} /2 . 

The theorem is thus clearly true when L 2p > 1 as soon as C p ^ 2, and this finishes 
the proof of (J3.3D . 

We now proceed to show that there exists C p (independent of M and k) such 
that (13.3[) holds, and decompose the proof of this fact into two steps. 



Step II. 1. We assume first that L 2p ^ 1, and consider again the martingale M as 
constructed in step 1.2. Since (M) 2 = 1, we know from step 1.1 that 

(3.9) d l (M 2 ,^)^C p Ll / p ( - 2p+1) . 

Let 

and observe that 

We have 

di(M 2 ,$) ^d 1 (M 2 ,M 2 ) + d 1 (M 2 ,« > ). 
The first term on the r.h.s. is smaller than S[|AM(r)|], and we have seen in (I3.7[) 

l/(2p) < 7-1/C 



M 2 = M T + S x _ (m)t _ , 
M 2 = M 2 -AM(r)l (M)l>1 . 



that this is smaller than L 2p p ' ^ L 2 p . Using also (I3.9|) , we obtain 

(3.10) di, fc (M 2 ,*)<di(F 2 ,*) < (C p + l)L 2 / (2p+1) . 

Let / € Cft(K,E) be such that ||/'||oo < 1 and ||/"||oo < fc. We will show that 

(3.11) \E[f(M 2 )} - E[f(M T )}\ < ~ (L^ + IKM)! - 1|U) . 
Indeed, since / € Cf(R,M) and ||/"||oo ^ k, we have 

|^ [/(M 2 ) - f(M T ) - (M 2 - M r )/'(M T )] | < ^ [(M 2 - M T ) 2 ] . 
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But 

E[(M 2 -M T )f'(M T )] = ^[ j B 1 _ (m>t _/'(M t )] 

= S[B[Bi_ Wt _ I J- t ]/'(M t )], 
and JS?[Bx-(m) - I ■?>] = since B and M are independent. On the other hand, 

E [(M 2 - M T f] = E^B^m)^ f] = E[l - (M) T -}, 
and we have seen in step 1.2, while treating the term 03, that 

E[1-(M) T -]<|KM) 1 -1|| 1 + L^. 
As a consequence, (|3. 1 1|) is proved, and thus 

(3.12) d 1)fc (M T ,M 2 ) < \ (l\I; + \\{M) X - lUx) . 
We now show that 

(3.13) di^Af^MO^IKAOi-lUi 
using the same technique. We write 

\E [/(Afi) - f(M T ) - {Mx - M T )f(M T )]\ < ^E [(Ah - M T f] , 

and observe that 

E [(Ah - M T )f'{M T )] = E[E[{M t - M T ) \ F T ] f(M T )] = 0, 

since M is a martingale and r a stopping time. On the other hand, we have seen 
while treating the term a\ in step 1.2 that 

E [(Ah - M T ) 2 ] <||(M)i-l||i, 

and thus ()3.13[) is proved. Combining (|3. 101) . Q3.12J) and (|3.13p . we thus obtain 

d ltk {Mi,*) < (C P + l) < 2p+1) + \h%* + k\\(M)i - l||i, 
and this proves (|3.3[) for L 2p ^ 1. 

S*iep //.,§. We now conclude by considering the case when L 2p > 1. We learn from 
step 1.3 that 

d ljfc (M 1 ,$)<2 + ||<M) 1 -l||l /2 . 
Since for any x ^ 0, we have yfx ^ 1 + a;/2, we thus obtain 

di, fc (Mi,*)<3 + i|KAf)i-l||i, 
and thus relation f|3.3[) holds when L 2p > 1, provided we choose C p ^ 3. □ 

4. The random walk in random conductances 

Let < a < (3 < +00, and tt = [a,(3] B . For any family u = (w e ) e£ B G £1, we 
consider the Markov process (X t ) t ^o whose jump rate between x and a neighbour y 
is given by w X]2/ . We write P" for the law of this process starting from x € Z d , E£ 
for its associated expectation. Its infinitesimal generator is L u defined in (II. ip . We 
assume that the (w e ) e6 B are themselves i.i.d. random variables under the measure 
P (with associated expectation E). We write P = PPg for the annealed measure. 
It was shown in [KV86] that under P and as e tends to 0, the process y/eX e -x t 
converges to a Brownian motion, whose covariance matrix we write A (in SS04 , 
it is shown that under our present assumption of uniform ellipticity, the invariance 
principle holds under Pq for almost every u>). 
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Let £ £ R d be a vector of unit L 2 norm. The purpose of this section is to give 



sharp estimates on the /c-Kantorovich distance between £ • X t /yt and &&(£), where 
we write $ CT to denote the cumulative distribution function of a Gaussian random 
variable with variance a 2 , and <r(£) = £ • At;. 

Theorem 4.1. For any S > 0, there exists a constant C (which may depend on 
the dimension) such that for any k ^ and any £ of unit norm, one has 



(4.1) 



di. 



V 



$ CT(0 ^c(fcvi)vr) 



for some q ^ 0, where in the l.h.s., £ • X t /V stands for the distribution of this 
random variable under the measure P, and where ^ q ,s was defined in (|1.3|) . 

Remark 4.2. When d ^ 3, the exponent of decay in (|4.ip can thus be made ar- 
bitrarily close to 1/2, and this is the exponent one gets when considering sums of 
i.i.d. random variables with finite third moment. 

Remark 4.3. By the same reasoning, one can also prove that there exist constants 
C (which may depend on the dimension) and q such that, for any £ of unit norm, 
one has 



di 



where again £ • X t 
measure P. 



V 



$ 



<r(£) 



<c 



r 



\0gl(t) t-V8 

logf (t) t-W 

£-1/4 



if d = l, 
if d = 2, 
if d = 3, 
if d > 4, 



£ stands for the distribution of this random variable under the 



Proof of Theorem \4-l\ The argument follows Moll] closely. We first treat the case 
d ^ 2. For (U > 0, we decompose £ ■ X t as M M (i) + R^t), where M^t) and -R M (£) 
are defined in [MoTTl (3.7)-(3.8)] respectively. Under the measure P (and for the 
natural filtration associated to X), M^ is a martingale with stationary increments, 
and we write 07, = E[M M (1) 2 ]. We have 



dijkl'^i,* 



v^ 



o-(0 



< di 1 



£•** M„(t) 



• di a 



M„(i) 



$ n 



■di, fc ($ 



°> >*"«)/ 



Vt ' Vt J ' V Vt 

with the understanding that random variables stand in place of their respective 
distributions under the measure P. Let us write the three terms in the r.h.s. above 
as 61 + 62 + 63, and proceed to evaluate each of these terms for the specific choice 
fi = 1/t. Considering (|2.5|) . we can bound the term b\ by 

E[(i? 1/t (t)) 2 ]V2 



di 



Vt ' Vt 



Ri/t(t)\ 



Vt 



< 



Vt 



and [Moll, Proposition 3.4] gives us that 



E[(R 1/t (t)) 2 } < C 



log 9 + (t) ifd = 2, 
1 if d > 3, 



for some constants C and g. 

To handle the term 63 , consider a standard Gaussian random variable A/". Then 
ctA/" has $0- as its cumulative ditribution function, hence 



di, fc (*„,$„>) < di (*„*•) < £[|<r.A/"- cr'AT|] = £[|JV|] |a - a'| 
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Since S[|A/"|] < 1, the term 63 is bounded by \a X j t — cr(£)|. It is shown in [GOlli 
Theorem 1] (and recalled in [Molll Proposition 3.3]) that 



Wi/t - cr(0| < C 



log^(i)^ 1 if d = 2, 

i- 3 / 2 if d = 3, 

log + (t) t- 2 if rf = 4, 

t- 2 if d 5s 5, 



which is much better than what we need for our purpose. 
We finally turn to the term b 2 . For any p > I, we introduce 



L 2p (t) = I E 



E l AM i/*WI 



2p 



O^s^t 



Theorem 13.21 tells us that if L 2p {t) ^ 1, then 62 is smaller than 

(M 1/t \- 1/2 



C p + -)(L 2p W) 1/(2p+1) + (fcVl) 



/t/t 



t 



°~i/t 



Replacing the exponent 4 by 2p leaves the proof of |MoII[ (3. II)] unchanged, and 
ensures that 

log q + (t) t-P +1 ifd=2, 

t- p+1 if d> 3, 



Lap(t) ^ C 



for some constants C and q depending on p. In particular, it is always true that 
L 2p (t) tends to as t tends to infinity. We fix p large enough so that 



(4.2) 



p-1 1 
2p+l > 2 



With such a choice for p, we have CL 2p (i)) 1/(2p+1) = o^ 1 / 2 ). 
On the other hand, |Molll (3.10)] ensures that 

logl^)*- 1 / 2 ifd = 2, 



(M x/ 



'/' 



fi/t 



^C* 



log+(t) r 



if d = 3, 
iid>4. 



Since 



(M- 



i/t/t 



&i/t 



€ 



(M; 



l/t)t 



V\/t 



this finishes the proof of Theorem 14.11 for d ^ 2 and i large enough, and it is easy 
to see that the l.h.s. of (|4.1[) is bounded for smaller t. The one-dimensional case is 
obtained in a similar way, following [Moll, Section 9]. □ 

5. HOMOGENIZATION 

We consider the discrete parabolic equation with random coefficients 

L^u on K + x Z d , 



(DPE W ) 



du 

~dt 



u(0, 



/ on Z d , 



where / : 7L d — >• R, L w is the operator defined in (|1.1|) . and by L u u(t,x), we 
understand L u u(t, -)(x). Note that L w is the discrete analog of a divergence form 
operator. 

For a fixed u £ ft, we say that u is a solution of (|DPE^p if it is continuous 
on [0, +00) x Z d , has continuous time derivative there (in other words, u(-,x) is in 
C 1 (R+,R) for every x G Z d ), and satisfies the identities dislayed in (|DPE"|) . 
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Proposition 5.1. For any lo £ f2 and any bounded initial condition f , there exists 
a unique bounded solution u of (IDPE"jl . and it is given by 

(5.1) u(t,x) = Klf(Xt)}. 

This is a very well known result. Checking that f|5. If) is indeed a solution is a 
direct consequence of the definition of the Markov chain. To see uniqueness, take 
u a bounded solution of (IDPE"j) . Letting M s — u(t — s,X s ), one can show that 
(M s )o^ s ^t is a martingale under P" for any x £ Z d , and as a consequence, 

u(t,x) = E£[M ] = E£[M t ] = E£[fi(0,X t )] = E£[/(X t )], 

which is the function defined in (15. ip . 

For a symmetric positive-definite matrix A, we consider the equation (jCPEft 
given in the introduction. We say that u is a solution of (jCPEp if it is continuous 
on R + x R d , has a continuous first derivative in the time variable and continuous 
first and second derivatives in the space variable on (0, +oo) x R d , and satisfies the 
identities dislayed in (jCPEp . 

Proposition 5.2. For any bounded continuous initial condition f , there exists a 
unique bounded solution u of (|CPE[) , and it is given by 

(5.2) u(t,x)=E x [f(B t )], 

where, under the measure P x , B t is a Brownian motion with covariance matrix A 
that starts at x. 



Again, this result is standard. It is proved in the same way as Proposition 15.11 
with the help of Ito's formula. 

Remark 5.3. The boundedness assumption in Propositions 15.11 and 15.21 could be 
changed for being subexponential. More precisely, let / : Z rf — > R be such that for 
any a > 0, \f(x)\ — 0(e a ^). Then there exists a unique solution u of (|DPE"[) such 
that, for any a > and any t ^ 0, sup s<t |u(s,a;)| = 0(e a ^). The boundedness 
condition was chosen merely for convenience. 

We now define rescaled solutions of the parabolic equation with random coeffi- 
cients. For a bounded continuous function / : R d — > R, we let vS e ' be the bounded 
solution of (|DPE"[) with initial condition the function x i-» f(ex), and for any t ^ 
and x € R d , we let 

(5.3) ue(t,x) =u^(e-% Le _1 *J) =Ef E - 1:c j[/(sX e - 2t )]. 

It is well understood (see for instance |BLP1 Chapter 3]) that the probabilistic 
approach yields pointwise convergence of u 6 to the solution of the homogenized 
problem. The following result is folklore (see also [LeOl] where the homogeniza- 
tion of random operators in continuous space is obtained using the probabilistic 
approach). 

Theorem 5.4. There exists a symmetric positive-definite matrix A (independent 
of f ) such that for every t ^ and x £ M. d , we have 

(5.4) u € (t,x) >u(t, x), 

where u is the bounded solution of (jCPEI) with initial condition f . 

Proof. Let us write (9 X ) to denote the translations on f2, acting according to 
(6 X w)j,, 2 = 0Jx+y.x+z- The distribution of X under P£ is the same as the one 
of X + x under P X " (both are Markov processes with the same initial condition 
and the same transition rates). Using this observation in (J5.3I) . we obtain that 

u £ (t,x)=-E e o lc - lxiU [f(eX £ -2 t +x £ )}, 



14 JEAN-CHRISTOPHE MOURRAT 

where x £ — e\e~ l x\. 

Since the measure P is invariant under translations, u e (t,x) has the same distri- 
bution as 

(5.5) E%[f(eX s - H +x s )]. 

It is proved in [KV 86, DFGW89] that for some symmetric positive-definite A (in- 
dependent of /), the quantity in (|5.5[) converges in probability to E [f(B t + x)] as 
e tends to 0, where B is a Brownian motion with covariance matrix A. □ 

Remark 5.5. It would be interesting to replace the convergence in probability in 
(I5.4[) by an almost sure convergence. Note that almost sure convergence for x — 
is equivalent to an almost sure central limit theorem for the random walk, and this 
is proved in [SS04 . Theorem 15.41 contrasts with for instance |JKQ1 Theorem 7.4], 
where weak convergence of an analogue of u £ is proved, but for almost every envi- 
ronment. 

We start the proof of Theorem 11.11 with two lemmas with a Fourier-analytic 
flavour. 

Lemma 5.6. Let Z be a random variable following the distribution v , M be a 
standard d-dimensional Gaussian random variable independent of Z , and a > 0. 
We have 



where 



and 



E[f(Z + aN)] = (2 7 r)- d ^exp ("^-) fiOHO d& 

fc) = Je**f(x)dx, 
i>(f) = J e -* x Av(x). 



Proof. Let us write 



Note first that 



9a{ - x) ~ (2^x 2 ) rf / 2 6XP V 2<r 2 



(5.6) g l/a (x) = exp \-L± J = (2^ 2 ) d ' 2 g a (x). 

The distribution of Z + <jN has a density (w.r. to Lebesgue measure) at point z 
which is given by 

g a (z - x) du(x) *P (27ra 2 y d/2 1 9i/a(z - x) dv(x) 



= (2na 2 )- d / 2 J e*<'-*>g 1/t ,(Z) d£ du(x) 

= (2iro*)- d ' 2 J e*'g 1/a (Z)u(t) <%. 

As a consequence (and using the fact that / and v are bounded), we have 
E[f(Z + o-N)} = (2^a 2 )- d ' 2 J f(z)e*>g 1/a (t)fi(Z) d£ dz 

= (2na 2 )- d ^ 2 J g 1/a (Of (OHO <%■ 
Since 

(2na 2 r d ^ 2 9l/ A0 = (27r)- d exp (-^) , 
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□ 



Lemma 5.7. For any integer m, there exists a constant C such that if if the weak 
derivatives of order m of f are in L 2 (M. d ), then 



{2*y d f(i + \e m ) |/(d 2 d^\\f\\i+j2\\d xT f\\i 

J 3=1 

Proof One has 

(-^) m / = %/ 
in the sense of distributions. By Parseval's theorem ( }RS[ Theorem IX.6]), since 
d a? f is assumed to be in L 2 (R d ), so is cX^f, and |!%/|| 2 = (2n) d / 2 \\d xT f\\ 2 . 
Hence, 



/' 



\^r\f(0\ 2 ^ = (2^r\\d xT m, 



and as a consequence, 



\e m \f(0\ 2 d^(2n) d J2P, T fW 



One also has ||/||2 = (27r) d/ ' 2 ||/||2, so the lemma is proved. □ 

Remark 5.8. In fact, as the proof reveals, there is a converse to the lemma: if 

(l+fc| am )l/(Ol a d£ 

is finite, then all the weak derivatives of / up to order m are in L 2 (M. d ). 
Proof of Theorem ] 1.1[ Let t > 0. We saw in the proof of Theorem 15.41 that 

E[u s (t,x)} = EEl^^lfieX^+x,)} 
= E[f(eX e - 2t +x e )], 

where in the last line, we used the fact that the measure P is translation invariant, 
and we recall that we write E for EEq and x £ for eL e_la; J • Note that 

d 

\E[f(eX e -2 t + x e )] - E[f{eX s - H + s)]| < ^ W d *J\\°° £ ' 

which is the first term in the r.h.s. of (|1.2p (a "lattice effect"). We now focus on 
studying the difference 

\E[f(eX E - H + x)}-E [f(B t + x)}\, 

where we recall that Eo[f(B t + x)] = Ei x [f(B t )] = u(t, x). Possibly replacing / by 
/( • + x), we may as well suppose that x = 0. Let a > be a small parameter, Af 
be a standard d-dimensional Gaussian random variable, independent of everything 
else, and write /( = f(Vt ■). Since ft is bounded and continuous, we have 



ft [ -^Xe-n + *M 



(5.7) E[f(eX e - H )]=E f t l— X e - H \ =ME 

Similarly, 

(5.8) Mf(Bt)] = Eo[/(V«Bi)] = E [/tCBi)] = lim E [/ t (Bi + aAf)] 



i() 
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where we slightly abuse notation by using the same N to denote a standard Gauss- 
ian (independent of everything else) under both the measures E and E. Let us 
write v e for the distribution of 

£ 

~7t 



-X e -2 t 



under the measure P, and vq for the distribution of B\ under Eo. Note that 

e £ ■ X e - H 



i> E (0=E 



exp i|£| 



V*l£l 



The function M. — >• R, x n- e l '^' x has first derivative bounded by |£| and second 
derivative bounded by |£| 2 .In view of (|2.4j) . we obtain from Theorem 14 . 1 1 that 



MO - Wl< c|£l (|e| V 1) V q ,s - 



Using Lemma 15.61 we thus obtain that 



E 



/ t ( ~r t Xe-H + oU 



<(2tt) 



-d 



exp 



a 2 |£| 2 



MO 1^(0 - A) (01 df 



€C* 



„.<5 



/*(0 Kl(|£|vi)d£, 



where C does not depend on a. Since /*(£) = t d ^ 2 f(£/Vi), we can perform a 
change of variables on the integral underbraced above, and bound it by 

(t + Vt) I f(0 (|£l 2 + i)d£. 



Let m = [d/2\ +3. By the Cauchy-Schwarz inequality, the integral above is bounded 
by 

1/2 , . ? \ !/ 2 



(l^ + l) 2 

2m ^ 



(i + iei 2m ) /(0 



cli 



1+ICI 2 

Since 2m — 4 > d, the first term of this product is finite, while Lemma 15 . 71 gives us 
that the second term is bounded by 

(STr^lll/ll^ + ^ll^/H 2 ) <(27r) d / 2 f ||/|| 2 + ^||%/|| 2 



We have thus proved that 



E 



f t [ ~r t Xs-H + ^ 



- EolMBi. + aM)] 



<c(t + Vi)[\\fh + J2\\d xT fh] * q , s 



Taking the limit a — > and recalling (J5.7I) and (|5.8|) . we obtain the announced 
result. □ 
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6. Heat kernel estimates 

The heat kernel pf (x, y) is defined so that (t, y) i— > pf(x, y) is the unique bounded 
solution to (|DPE"[) with initial condition / = l x . The heat kernel is symmet- 
ric: pf(x,y) = p%(y,x), and by translation invariance of the random coefficients, 
EW(x,y)]=E\fi(p,y-x)]. 

The aim of this section is to prove Theorem 1 1.21 In order to do so, we will need 
a regularity result on the averaged heat kernel. For / : Z d — > R and 1 ^ % ^ d, we 
write 

Vif(x) = f(x + ei)-f(x), 
where (ej)i^i^ is the canonical basis of M. d . The following result was proved in 
[CNOOal Theorem 1.4], and then elegantly rederived in JDD051 (1.4)]. 

Theorem 6.1 ( [CNOOal IDDOB] ) ■ Let 

(6.1) q t (x)=E\p?(Q,x)]. 

There exist C,c\ > such that for any t > and any x £ Z d , one has 

_ C ( f\x\ 2 

\VM*)\ < ^piTT^ exp \~ Cl \t A N 

We also recall the following upper bound on the heat kernel, taken from [De99, 
Proposition 3.4] (see also [CKS87I Section 3] for earlier results in this context). 

Theorem 6.2 ( Dc99 ). (1) There exist constants C,c such that for any t ^ and 

any x E Z d . 

C 



#(M<j-7^ exp (-AstO*)). 



where 



D,[.r) = .r|ar,inll ( ^1 + t I ^ 1 + ^- - I 



(2) In particular, there exists C2 > such that for any x E Z d 
Pt(0,x)^ I ^c^f-cJ l -^fA\x\ 



Proof of Theorem ] 1.21 We decompose the proof into three steps. 

Step 1. Possibly lowering the value of C2 > 0, we have that for any x E M. d , 

(6.2) p l (0,x)^Cexp(-c 2 \x\ 2 ), 

(6.3) ^°' ° (a) < Cexp (-c 2 |x| 2 ) (1 < t < d). 

Equation (|6.2p and part (2) of Theorem 16.21 thus ensure that (possibly enlarging 
C), 

(6.4) \e~ d g e - a (Le _1 xJ)-3J 1 (0 > !t)| < Cexp (-c 2 (|:r| 2 A |e _1 a;|)) . 

Moreover, Theorem 16.11 remains true if we lower the value of the constant c\ > in 
such a way that c 2 ^ c±/2\/d. 

Step 2. We now show that there exist c > (independent of 5) , es > and Ca such 

that, for any e ^ e^ and any iGl 1 ', one has 

(6.5) 



\e- d 



q e -*{\s- l x\) -^(0,x)| ^ C, (*^(e 2 )) 1/(d+3) exp (-c(|x| 2 A le" 1 *])) 



Let / be a positive smooth function on R d with support in [— 1, l] d and such that 
J f = 1. We define, for any r > 0, the function f r :x>-^ r~ d f(r~~ 1 x). 



l,x 
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Let u^ be the bounded solution of (|DPE"P with initial condition f r (e •) (we 
keep the dependence of u^ £ > in r implicit in the notation). By linearity, we have 



,(<0 



{t,x)= ^2Mez)tf(z,x). 



Letting u e (t,x) = u^(e 2 t, [e x x\), we obtain 

(6.6) u e (t,x) = Y^ fr{ez) p"- H (z, [e~ x x\] 



Let u be the bounded solution of (jCPEp with initial condition /,.. Observing the 
proof of Theorem ll.il we get that for any S > 0, there exists C such that 

(6.7) |E[u s (l,s)]-«(l,a:)| 



d „ 

< E H^Mloo e + C * q Ae 2 ) I \f r (0\ (\e + 1) d£ 



Scaling relations ensures that ||9x,/r||oo is bounded, up to a constant, by r ( d+1 ), 
while f r (0 — /(*"£)• As a consequence, 



/, 



UO ie| 2 d^-r 



„-d 



-(d+2) 



f(0 ICI 2 ^, 



and the integrals on the r.h.s. are finite since / is smooth (see Lemma [577]) . To sum 
up, for some constant C and any r ^ 1, we have 

(6.8) |E[u £ (l, a:)] - «(1, a;)| ^ C (e r- {d+ ^ + tf, |4 (e 2 ) r~ (d+2 ^ . 

The solution u can be represented in terms of the heat kernel as 

u(l,x) = I f r (z)p x {z,x) dz = p 1 (0,x) + / f r {z)(j) x {z,x)-p x {0,x)) dz, 



where we used the fact that J/ r = 1. For z £ M d such that ||z||oo ^ T ^ 1 and up 
to a constant, \p 1 (z,x) — p 1 (0,x)\ is bounded by re~ C2 ' x ' by (|6.3p . Since f r has 
support in [— r, r] d , we arrive at 

(6.9) |u(l,a:)-Pi(0,a;)| < C rexp (-c 2 |x| 2 ) . 
On the other hand, if z £ Z d is such that ||2||oo ^ s~ l r, then 

|E[p£_ 2 (z, Le _1 a;J)] -q e -2(Le _1 a;J)| < de' l r sup |V<g s - 2 ([e -1 a;J -*)| 

||*||oo<£~ 1 »- 

We now argue that there exists C3 > (independent of 5) such that, uniformly over 
r ^ 1 and x € R d , one has 

(6.10) sup IViJfe-^Le^xJ-^K-jfieo^f-CBdajI'Ale- 1 *!)]. 

I|z||oo<£ _1 r £ 

l<i<d 



Theorem 16.11 tells us indeed that the l.h.s. of (I6.10[) is smaller than 

ILrtj-zl 2 



c 

^TT ex P 



— C\ inf 

l<i<d 



A \[e~ L x\ -z\ 
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For any r ^ 1 and ||x||oo ^ 2, the infimum above is larger than 

\x\ 2 A le" 1 ^! 

so (|6.10p holds in this case, with c 3 = c\j2\fd. To control smaller values of ||x||oo, 
it suffices to enlarge the constant C in f|6.10[) . To sum up, we have shown that 

|E|j£_ a (*, Le" 1 ^)] - q e -2([e- l x})\ ^ C e d rexp [-c 3 (M 2 A le -1 ^)] . 

In the sum on the r.h.s. of (I6.6[) . only C(e~ 1 r) terms are non-zero, and ||/||oo ^ 
r~ d , so 



E[u e (l,x)] ~ J2 /r(«) ff E -»(Le _1 a:J) 



«e2 



< C rexp [-C3 (|x| 2 A |e _1 x|' 



Observe also that 

This is a Riemann approximation of J f = 1, hence 



£ d £/ r ( £ z)-l 

zSZ d 



r 



and we are thus led to 

(6.11) |E[u s (l,x)] - e~ d g e - 2 (Le- 1 xJ)| < C (rexp [-c 3 (|x| 2 A |e- x a:|)] + J) . 

Combining (|6.8[) . (16.91) . (|6.1ip and the fact that C2 ^ c 3 = c\/2y/d, we obtain that 
up to a constant, 

|e" d fe-»(|e"Vl)-Pi(O,30| 
is bounded by 



r d+i 



r d+2 



rexp [-c 3 (|x| 2 A le" 1 ^!)] + - 



uniformly over r sj 1. Since for e small enough, one has e ^ ^q,s (=); the above is 
bounded, up to a constant, by 



(6.12) 



%^P+ rexp [-c 3 (M 2 A |rt|)] , 



uniformly over r ^ 1. Choosing 

r d+3 = *^( £ 2 ) exp [c 3 (|x| A \e- l x\ A M £ )] , 

where 

log(*,, 5 (e 2 )) 



M F 



c-.i 



is here to ensure that r ^ 1, we obtain that the expression in (]6.12[) is smaller than 

1 



(^( £ 2 )) 1/(rf+3) exp 



c 3 1 



rf + 3 



(|x| A |e _1 a;| AAf E ) 
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This proves (|6.5p when |x| A |e _1 2:| ^ M e . Otherwise, we use the bound (|6.4p , 
together with the fact that C2 ^ C3, to get 

I e-^g^de -1 ^)- Pi M | 

^Cexp(-c 3 (|x| 2 A|e^ 1 a;|)) 

^ e xp(- C 3(l-^)(|^A|rt|)--^M £ 

,C C (* 9 ,,(e 2 )) 1/(d+3) exp (-c 3 (l - -±-^ (N 2 A | e -i*|)) . 

Hence, (|6.5p holds also in this case, and we can always choose c = 03(1 — l/(d + 3)). 

Step 3. We now extend the result to any time t > 0. The heat kernel of the 
continuous operator satisfies the scaling relation 

p t (o,x) = t- d/2 p 1 (o,x/y/i), 

while we can write 

e- d fc-atde- 1 *]) = * _d/2 (eM)" d ^/^-(L^M)- 1 (*M)J)- 
For es and C,5 given by step 2, as soon as e/y/i ^ ej, one thus has 



C S ( / £ 2XxW+3) 



< ^ I **■* It ' ' exp 



- ( (M! A | £ -^i 



which is the claim of the theorem. D 

7. HOMOGENIZATION OF ELLIPTIC EQUATIONS 

In this last section, we state and prove the counterparts of Theorems II. II and II. 21 
for the homogenization of elliptic equations. For / : M. d — > M bounded continuous, 
we consider the unique bounded solution of 

(DEE^) (e 2 - L u )v {e) = e 2 /(e •) on Z d . 

Using integration by parts, one can check that 

p+ao 

(7.1) v {e \x) = e~* u {e) {e-H,x) di, 

Jo 



where vf- e > is solution of (|DPE"[). For x £ K d , we let v E (x) = v^([e 1 x\), so that 



r+00 
(7.2) v E (x) = / e~* u e (t,x) dt 

Jo 

The function v e converges pointwise, as e tends to 0, to v the bounded solution of 
(CEE) (l--V-AVJv = / onK d , 

and one has 

-t 



(7.3) v{x) = / e _t u(i,x) d£, 

Jo 

where u is the solution of ([CPEp . Equipped with the representations (|7.2p - (|7.3p . it 
is straightforward to derive the following result from Theorem ll.il 
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Theorem 7.1. Let m = [_d/2\ + 3 and 6 > 0. There exist constants C$ (which 
may depend on the dimension) and q such that, if the weak derivatives of order m 
of f are in L 2 (M. d ), then for any e > and x £ M. d , one has 



\E[v e (x)] - v(x)\ < Yl H^i/IU e + C s [\\fh + J2 W&Tfh *<,,* (e 2 ) • 

3=1 \ 3 = 1 j 

Remark 7.2. Note that on the other hand, it does not look so simple to deduce 
Theorem 1 1 . 1 1 from Theorem 1 7. II A possibility for doing so may be to try to devise 
a quantitative version of [Kal Theorem IX. 2. 16]. 



One can also consider the Green function G£(x, y), the unique bounded function 
such that 

(e 2 -L")G^(x,-) = l x . 

Letting G(x, y) be the Green function associated to equation (|CEE[) . we can write 
the counterpart of Theorem 11.21 

Theorem 7.3. Let d Js 2 and 5 > 0. There exist constants c > (independent of 
8), q, Cs such that for any e > and any x € eL d \ {0}, one has 



(7.4) \e 2 - d E [G£(0, e^x)] - G(0, x) 

Cs 



€ 



Id— 2 



* 



q ' 5 l 'd 2 



l/(d+3) 



- c \' x \ 4. e _c ' £ la: ' 



When d = 1, there exist C, c > such that for any e > and an?/ a; £ eZ, one /ias 
leEfG^O.e- 1 ^)] -G(0,x)\ < C* 



g l/8 e -c|x| +e -c|e ^1 



Remark 7.4. The orders of magnitude, as e tends to 0, of the r.h.s. of (|7.1[) and 
(|7.4|) . are given respectively by (|1.5p and (|1.6|) . 

Proof. Our starting point is the fact that 



/■+oo 

G?(x,y)=e- 2 e-t p u s - H (x,y) dt, 

Jo 



while 



G(z,y) 



+00 



e * PtfolO dt 



Recall first that Theorem 11.21 ensures that there exist c > 0, Cs,es > such that 
whenever t ^ (e/es) 2 , one has 



(7.5) \e- d E^. H (0,[e- 1 x\)]-p t (0,x)\ 



The difference of interest 



,2-d 



S rf / 2 "^ \t ,exp 



E[G£(0, |e _1 ajj)] -G(0,x) 



A |£ _1 a;| 



is bounded by 
(7.6) 



+ ■00 



e"* Ie- d 



I[pgL, t (0,Le- 1 a!j)]-p t (0,»)| di. 



Let 77 = (e/e^) 2 V (e|ar|). If t ^ ry, then the integrand above is bounded, up to a 
constant, by 



e ,i/(d+3)/e 



^r 1 r 

£d/2 9,5 



f 



exp 



21 



/ 
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In order to control the integral in (|7.6[) , it thus suffices to bound the following three 
quantities: 

12- 






Fr 

-c 



elf. 



(7.8) r e -' J E[p£L at (0,Le- 1 xJ)]d*, 



(7.9) P p t (0,x) dt. 

Jo 



We start with the integral in (|7.7p . which is the only non- negligible one. To begin 
with, note that for any 7, a change of variables gives us the identity 

(7.10) j o "—e-^ 2 /* dt = \x\^-< j Q — r e~ c ' s ds, 

and moreover, provided 7 > 1, 

+00 -slzl 2 /.l/|x| -c/2s /-+oo _ c / 



- c / s d S < e -VA/2 1 da + e-W / ds 



« 7 Jo s7 A/|x| s7 

(7.11) < Ce- cN/2 , 

for some large enough C (and c ^ 2). We have thus shown that, for 7 > 1, 
/•+00 —t 

(7.12) / ^e"^ '/* dt < Clx^-^e- ^!/ 2 . 

When <i ^ 3, we have \E' ? ^(m) = u 1 / 2 "* 5 , so that the integral in (J7.7I) is bounded, up 
to a constant, by 

(7.13) |*f-' *#*«> (_fl) e -cW/2. 

When d = 2, the argument requires some minor modifications, due to presence of 
a logarithmic factor in ^5,5. One should consider instead integrals of the form 

i^^ogi (t/e 2 ) e-^ 2 /t dt = M 2 - 2 - / + °° ^log 9 ! ( S |x| 2 / £ 2 ) e"^ ds, 
Jo l ' Jo s 

for some </ ^ and 7 > 1 (in fact, 7 = 1 + 1/20). This last integral is bounded by 



L 



£ 2 /M 2 e -s\x\ 2 /-+OC e -s\x\ 



s7 ./e 2 /M 2 s7 



/s ds + / log 9 (s\x\ 2 /e 2 ) e- c l s dt 



For the first integral, (|7.11[) gives us an upper bound. Inequality (|7.1ip also enables 
us to bound the second integral, using the fact that 

log 9 ' (s\x\ 2 /e 2 ) < 2 9 ' (log 9 ' (\x\ 2 /e 2 ) +log 9 ' («)) . 

These observations thus guarantee that (I7.7[) is also bounded by (17.131) when d = 2. 
We now turn to the evaluation of the integral in (|7.8p . Since, for z Js 0, one has 
arsinh(z) = log(z + y/1 + z 2 ) ^ log(l + z), and using part (1) of Theorem 16. 2[ one 
can bound the integral in (|7.8I) (up to a constant) by 

/Vc XP ( H rt|i„ g (i + Jg|i))d ( . 

A change of variables shows that this is equal to 

(7.14) ^ £ -d J exp (_| e -i x | log (i + 1/s )) dSj 

c Jo 
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where 

. cri ce7 

e\x\ \s L x\ 
Since we consider only x € eZ d \ {0}, the parameter rj is uniformly bounded, 
independently of the value of x and e. The integral in (I7.14[) is thus bounded (up 
to a constant) by 

e^larKl + ir 1 )" 1 *" 1 * 1 = \x\ 2 - d \e' 1 x^' 1 {1 + -q- 1 )-^' 1 ^ 

^C\x\ 2 - d exp(~c\e- 1 x\). 



This finishes the analysis of the integral in (|7.8[) , and there remains only to consider 
the integral in (I7.9p . This integral is bounded by a constant times 

rv t -*/* e -WH d* 

/o 

for some small enough c > 0. A change of variables enables one to rewrite this 
integral as 

(7.15) 

|.z;| 2 - d / u- d / 2 e- c ^ u du < |x| 2 - d exp (__£_) / u - d / 2 e- c / 2 " du. 

^r 2 = r^v-i-<-^-<c' 



'o 

Moreover, 



for some large enough C", uniformly over e > and a; <G eL d \ {0}. The r.h.s. of 
(|7.15[) is thus bounded by 

|x| 2 - d ex P (-^i) / C U - d/2 e" c / 2 " du. 

We thus obtained the required bound on (I7.9P , and this finishes the proof of Theo- 
rem O for d^2. 

For the one-dimensional case, the analysis must be slightly adapted. We need to 
bound the integrals appearing in (|7.7p , (I7.8[) and (17.91) . The analysis of the integrals 
in (J7.8I) and (I7.9[) can be kept without change, except that only the case x S eZ\{0} 
was considered above, while here we want to consider also x — 0. But this is a very 
easy case, since the upper bound t^ 1 / 2 on the heat kernels is integrable close to 0. 
As for the integral in (|7.7[) . it is equal to 

"+oo _* 



-1/ 



o * 7 



where 7 = 1/2 + 1/16 < 1. The integral above is uniformly bounded over x such 
that |jb| ^ 1. Otherwise, as noted in (J7.10I) . we have 

+00 -t , f+00 „-s\x\ 2 

e ~ c l-l 2 /* dt - M 2 - 2 ^ / e~ c ' s ds, 



V h ^ 



and we can bound the last integral by 

-i/M e -« ,, f+°° e - s / 2 



<? 



W/ !_d S + e-N/ 2 / ds, 



sl J\/\x\ sl 



where in the second part, we used the fact that for \x\ ^ 1 and s > |a; | x , we have 
s|x| 2 ^ |*|/2 + s/2. We have thus shown that the integral in (J7.7I) is bounded, up 
to a constant, by 

E 1 ^ (1^2-27 + 1 ) e -c|x| ) 
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uniformly over i£l, and this finishes the proof for d = 1. □ 
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